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SUMMARY 


A study has been made of the effects of varying the shape, solidity, 
and heat-transfer coefficient of thin wings with regard to their influ- 
ence on the torsional- stiffness reduction induced by aerodynamic heating. 
The variations in airfoil shape include blunting, flattening, and com- 
bined blunting and flattening of a solid wing of symmetrical double-wedge 
cross section. Hollow double-wedge wings of constant skin thickness with 
and without internal webs also are considered. The effects of heat- 
transfer coefficients appropriate for laminar and turbulent ilo w Qxe 
investigated in addition to a step transition along the chord from a 
lower to a higher constant value of heat- transfer coefficient. From 
the results given it is concluded that the flattening of a solid double 
wedge decreases the reduction in torsional stiffness while slight degrees 
of blunting increase the loss . The influence of chordwise variations in 
heat- transfer coefficient due to turbulent and laminar boimdary-layer 
flow on the torsional stiffness of solid wings is negligible. e e 
of a step transition in heat- transfer coefficient along the chord of a 
solid wing can, however, become appreciable. The torsional-stiffness 
reduction of multiweb and hollow double-wedge wings is substantially less 
than that calculated for a solid wing subjected to the same heating 
conditions . 


INTRODUCTION 


The influence of aerodynamic heating on the effective torsional 
stiffness of thin wings was investigated in reference 1. From the 
results presented, the loss of torsional stiffness in thin wings, rapidly 
accelerated to supersonic speeds, was found to be significant. Th s con- 
cision was based on calculations for a solid airfoil of double- wedge cross 
section with a constant heat-transfer coefficient. A hollow double-wedge 
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wing of constant skin thickness with a variable heat-transfer coefficient 
was also investigated. The magnitude of the loss in torsional stiffness 
was found to depend primarily on the chordwise mass variation of the air- 
foil cross section. 

In the present paper the influence on torsional- stiffness reduction 
of changes in the mass and temperature distributions of the airfoil cross 
section is studied. The changes in the mass distribution were effected 
in two ways: By considering various modifications of a solid double- 

wedge cross section, and by considering an airfoil of multiweb construc- 
tion. The changes in the temperature distribution (other than those 
which accompanied the change in the mass distribution) were effected by 
varying the heat- transfer coefficients. The effects of heat-transfer 
coefficients appropriate for laminar and turbulent flow were investigated 
in addition to a step transition along the chord from a lower to a higher 
constant heat-transfer coefficient. Details of the analyses are given 
in the appendixes. 


SYMBOLS 


A cross-sectional area 

a semiwidth of flat portion of modified wedge 

B arbitrary constant 

b modified wedge dimension (see fig- l) 

by web spacing 

c chord 

c m specific heat of wing material 

d distance from leading edge of wirg 

E Young's modulus 

G shear modulus 

h heat-transfer coefficient 

I moment of inertia 

J torsional-stiffness constant 


K 


constant 
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L = a + b 
M 


N- 


Pr 


R 

T 

Zff 


aw 


x 

a 

P 

7 

e 

n r 

e 

A 


thermal conductivity of air 

c. 

2 

Mach number 
Prandtl number 


Reynolds number 

average temperature through thickness of wing 

difference between initial and final adiabatic wall 
temperature 

wing thickness 

distance from midchord of cross section 

coefficient of thermal expansion 

dummy variable 

ratio of specific heats 

spanwise strain 

temperature-recovery factor 


T - T 

temperature parameter, — — — - 

^aw 


time parameter 


I = 


2d 

c 


p mass density of wing material 

a normal stress in spanwise direction 

T time 

$ torsional-stiffness-reduction parameter (see eq. (12)) 
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Subscripts : 


average 


adiabatic wall condition 


at midchord 


effective 


final 


laminar 


leading edge 


maximum 


initial flight condition 


total 


trailing edge 


turbulent 


free stream 


BASIC THEORY 


The torsional behavior of thin wings cun be strongly influenced by 
the presence of axial stresses. An expression for the effective tor- 
sional stiffness GJ e ff of a thin wing of arbitrary cross section 

subjected to axial stresses has been derived in reference 1 and is 


J eff _ j + / crr-dA 
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where 

GJ Q initial torsional stiffness 

a spanwise normal stress (tension positive) 

r distance from axis of twist to point in question 

When the wing is accelerated to supersonic speeds aerodynamic heating 
induces spanwise thermal stresses in the wing structure which in turn 
cause a change in the torsional stiffness as shown in equation (l) . The 
seriousness of this change depends on a variety of factors such as flight 
condition, heat-transfer coefficient, mass distribution, and the thermal 
properties of the wing material. 


Temperature Equations 

In order to determine the chordwise temperature distribution (and 
consequently the spanwise stress) in a wing accelerated to supersonic 
speeds, the following assumptions are made: 

The wing is assumed to be instantaneously accelerated from an 
initial Mach n umb er Mq to a final Mach number Mf at a constant alti- 
tude. This assumption is based on the conclusions reached in reference 1 
that the max imum loss of torsional stiffness for moderately rapid accel- 
erations is only slightly less than that calculated for instantaneous 
acceleration. The local free- stream temperature and Mach number are 
taken to be those of the free stream ahead of the wing. Thus, the free- 
stream temperature is constant throughout the flight and the heat- transfer 
coefficient h is not a function of time. The heat-transfer coefficient 
is considered to be symmetric about the chordline corresponding to a 
0° angle of attack. 

The initial wing temperature at time zero is assumed to be uniform 
and equal to the adiabatic wall temperature at the initial free- stream 
Mach number Mq. The adiabatic wall temperature is given by 


aw 


= Tjl + 


1 


T l r ( 7 - 1)M2 


( 2 ) 


where 

local free-stream. temperature 
r\ v temperature-recovery factor 



6 


Chordwise and spanwise heat conduction, as well as radiation effects, 
are neglected. Thus, a heat balance equation governing the chordwise tem- 
perature distribution can be expressed in terms of the average temperature 
through the thickness as 


h ( T aw - T ) = P c m | ^ (5) 

where 

T average temperature through thickness of wing 

T time 

p mass density of wing material 

c m specific heat of wing material 

t wing thickness 

Consequently, the average temperature at any point along the chord can 
be expressed in nondimens ional form as a function of a time parameter 

A and a ratio of wing thickness to thickness at the midchord — — as 

^max 

9 = T " T o = i _ g^/^raax (4) 

^aw 


where 

To 

*Tav 


and 


initial wing temperature 

difference between initial and final adiabatic wall 
temperature 


A = 


2hr 

pc m t max 


(5) 


From equation (5) it can be seen that for constant values of p, 

c m , and t_ 0 _, the chordwise variation of \ is proportional to the 
ui 7 max 

chordwise variation of h, which in turn is dependent mainly on the 
condition of the fluid flow over the airfoil surface. 



7 


Heat-Transfer Coefficients 


In order to study the influence of a variation in heat- transfer 
coefficient along the chordlength on the stiffness reduction, heat- 
transfer coefficients appropriate for turbulent and laminar boundary- 
layer flow need to be considered. For this analysis, flat-plate heat- 
transfer coefficients are assumed, and the Mach number on all surfaces 
is assumed to be the same as that of the free stream ahead of the wing 


For the case of turbulent boundary- layer flow, the heat- transfer 
coefficient at a point can be expressed by the empirical equation given 
in reference 2 as 


Hurb _ 


0 . 029 k a R d ^^Np r 1 / ? 


(6) 


where d is the distance from the leading edge to the point in question. 
By using the known identities for Mach number and Reynolds number, equa- 
tion (6) can be rewritten in terms of and may be expressed in 

terms of the heat-transfer coefficient at the midchord as 


^turb “ 


N;urb,c/2 

(f ) i/5 


( 7 ) 


Similarly, if the flow is laminar the heat-transfer coefficient can be 
obtained by the use of the empirical equation given in reference 2 and 
is expressed as 


1 lam ~ 




( 8 ) 


Equation (8) can also be rewritten in terms of the heat- transfer coeffi- 
cient at the midchord as 


hlam, c/2 

¥)‘ ,! 


( 9 ) 


Thermal-Stress Equations 

Once the heat- transfer coefficient is established and the tempera- 
ture distribution is known, the thermal stresses can be calculated on 
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the basis of elementary beam theory in which it is assumed that plane 
sections remain plane with no lateral restraint. The material properties 
are considered to be invariant with temperature and all deformations are 
assumed to be elastic . 

Since the doubly symmetric sections considered in this paper can 
undergo only a chordwise rotation and translation of the cross section, 
the spanwise strain can be written as 


e = K + Bx = | + a(T - T q ) 
where K and B are constants . 


( 10 ) 


Application of the equilibrium conditions of zero thrust and zero 
moment permits the evaluation of K and B in equation (10) . The 
resulting expression for the spanwise stress can be written as 


q 

EaZffaw 
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( 11 ) 


where 


■'av 



8 dA 


When the spanwise stresses on a free -end beam are due to temperature 
they are necessarily self-equilibrating and consequently the magnitude of 


the thermal- stress torque given in equation (l) 



is independ- 


ent of the axis of twist. For the doubly symmetric sections considered 
in this paper, however, the axis of twist s at the midchord of the 
cross section and for thin wings r may be approximated by x. 


By substituting equations (4) and (lli into equation (l) and per- 
forming the necessary integrations, the change in torsional stiffness 
can be expressed as 


^eff _ ± 
GJ 0 


Ea 

G 


^•av 



<t> 


( 12 ) 
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where the parameter $ is a function of the time parameter A and the 
shape of the cross section, hut it is not a function of the thickness 

ratio 

c 

lW ]QY 

Since the thickness ratio — - — , material properties, and change in 

adiabatic wall temperature ££T av are constant for an airfoil at a given 

flight condition, the percent reduction in torsional stiffness will depend 
solely on the magnitude of 4>. Furthermore, values of 0 may be compared 
directly for each of the examples to be considered in the following sec- 
tions of the paper. The influence of varying the shape and solidity of 
the cross section and of the distribution of the heat-transfer coefficient 
on the percent reduction in torsional stiffness will therefore be deter- 
mined by investigating the change in <P. 

AMLYSIS 


Effect of Variations in Cross-Sectional Shape 

In order to study the effects of varying the airfoil shape on the 
reduction of torsional stiffness in solid cross sections, the modified 
double-wedge wing shown in figure 1 was analyzed. The analysis was 
based on the assumption that the heat-transfer coefficient remained con- 
stant across the chord, and therefore A in equation ( 5 ) is a function 
of time alone . 

The temperature distribution and the thermal stress distribution 
are symmetric about the midchord and only a translation of the cross 

section can occur (i.e., B = 0 in equation (10) and / 0x dA = 0 

^A 

in equation (ll)). Expressions for the temperature and stress distri- 
butions over the cross section and the percent of torsional- stiffness 
reduction for the modified double-wedge wing are derived in appendix A. 


In order to show the quantitative effects of flattening and blunting 


of a cross section on the torsional-stiffness reduction, the parameter $ 

was evaluated for different degrees of flattening — and blunting ^~ e 

c t rnax 


over short time intervals. The results are shown in figures 2 to 5* 


The influence of the flattening only of a solid double-wedge cross 
- 0 on the torsional-stiffness reduction is shown in 

^max 


section 
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figure 2. The stiffness parameter is plotted as a function of the time 
parameter A for various degrees of flattening, ranging from the unmodi- 
fied double-wedge ^ = 0^ to the flat plate = ij. The curve for 

a. = 0 is identical to that given in reference 1. In the limiting case 
of — = i the wedge becomes a flat plate and, as would be expected, the 

c 2 

reduction in torsional stiffness is zero. It is evident from figure 2 
that flattening of a solid double wedge reduces the loss of torsional 
stiffness caused by aerodynamic heating. The greatest loss in torsional 
stiffness = 0.04-56) occurred for the unmodified double wedge. As 
shown in reference 1 this value of would correspond to a 50-percent 
reduction in torsional stiffness for a 5-percent-thick stainless- steel 
wing accelerated instantaneously from a Mach number of 0 to a Mach num- 
ber of 4 at an altitude of 50,000 feet. 


The influence of the blunting of a solid double-wedge cross sec- 
tion on the torsional-stiffness reduction can be studied by considering 
the flattened portion of the modified wing s./c to be equal to zero. 

In figure 5, $ is plotted as a function of A for various degrees of 

blunting. The degree of blunting is given bj the ratio of tip thickness 


t]_ e to midchord thickness t maY . 


When 

^max 


0, the reduction in tor- 


sional stiffness Is identical to that calculated for a solid double wedge 

( 5: = 0 in fig. 2). As approaches unity the double wedge becomes 

' c > ^max 

a flat plate and the torsional- stiffness-reduction parameter $ approaches 
zero. For a slight degree of blunting, however, there is an increase in 
the max imum value of 4>. In order to obtain a clearer picture of this 
variation, the maximum values of $ were plotted as a function of the 

blunting ratio in figure 4. It is evident from figure 4 that for 

4max 


slight degrees of blunting, the maximum torsional-stiffness reduction 
is greater than that obtained using the unb tinted double wedge. This 
increase in stiffness loss reaches a maximum of 11 percent for a blunting 
ratio of approximately 0.10. Blunting ratios greater than 0.20 are 
needed to decrease the maximum torsional-stiffness loss below that cal- 
culated for the unblunted wedge . 


In figure 5 the effect of combined blunting and flattening of a 
solid double-wedge section is shown. The stiffness-reduction parameter 
<D is plotted as a function of A for various degrees of blunting 

of the leading and trailing edges of a flattened wedge where 

^max 

a/c = l/4. The maximum reduction in torsional stiffness as a function 
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of the blunting ratio — is plotted in figure 4. In this figure the 

' t 'max 

max imum stiffness reduction for the blunted and flattened double wedge 
can be seen to be similar to that of the blunted wedge. The effect of 
flattening of the midchord, however, reduces the magnitude of the maxi- 
mum stiffness reduction by approximately 50 percent in comparison with 
the blunted wedge . 


Effect of Variations in Solidity 

In order to investigate the effects of variations in solidity on 
the reduction of torsional stiffness, a wing of multiweb construction 
such as that shown in figure 6 was considered. The results of an analy- 
sis of the torsional-stiffness reduction due to aerodynamic heating of 
such a multiweb wing were presented in reference 5 for a constant heat- 
transfer coefficient. The analysis, however, was based on the assump- 
tion that the temperature of the webs remained unchanged while the 
temperature of the skins asymptotically approached the adiabatic wall 
temperature. Thus, the temperature gradient and consequently the ther- 
mal stresses of the cross section continue to increase to an asymptotic 
value since the alleviating effect of heat conduction between skin and 
webs is not present. 

In the present analysis, the heat-transfer coefficient was con- 
sidered constant, but in order to allow the temperature of the skin and 
web material eventually to equalize, the webs were allowed to absorb 
heat. The temperature of each skin was considered to vary as that of a 
flat plate and can be written as 


9 = 1 - e 2ts / tmax (13) 


The temperature of the web material was considered to vary as if 
the web mass completely filled the multiweb and formed a solid double 
wedge. The temperature of individual webs is then established as the 
temperature of different chordwise locations on a solid double wedge so 


that 


-A 

e„ = i - e t/tmax 


(iM 


The average temperature in a typical bay of the multiweb can be 
expressed as 

2byt s B s + 


0 


^b^tg + ttntf 


(15) 
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In order to evaluate the spanwise stress at any point along the 
chord the total area of the web material is considered to vary contin- 
uously along the chord in proportion to the wing thickness. The expres- 
sion for 0 in equation ( 15 ) can then be substituted into equation (ll) 
for spanwise stress. Because of the symmetry of the temperature distri- 


bution about the midchord, 


0x dA in equation (ll) is equal to zero. 


The expression for the torsional- stiffness-reduction parameter $ 
of the multiweb is obtained in appendix B. In this analysis the tor- 
sional constant J Q was calculated on the basis that the webs were rigid 
in shear . 


In figure 7 ® is plotted as a function of A for various ratios 
of skin thickness to maximum chord depth and ratios of web area to skin 
area A w /A s . For a specific ratio of A^/Ag the maximum stiffness 

t K 

reduction increases with a decrease in the ratio of — — , whereas for 

tmax 

t Q 

a specific ratio of — — , the maximum stiffness reduction decreases with 

tmax 

a decrease in A^/Ag. When t w approaches zero the multiweb reduces to 

a hollow double wedge and the drop in torsional stiffness becomes zero. 
Note that the maximum value of 0 is somewhat lower than that calculated 
for a solid double wedge . 

In figure 8 a comparison between the torsional-stiffness reduction 
obtained from the present analysis and that obtained from the analysis 
based on the assumption of no heat conduction (ref. 3) into the webs is 
made for two multiweb wings. The wings are assumed to be instantaneously 
accelerated from an initial Mach number of 0.75 to a final Mach number 

of t.O at an altitude of 50,000 feet. For bo :h wings ^ ax and the 

web area to cover area A w /A s are taken to be, respectively, 0.03 

and 0.35* The skin thickness t G and ratio of web thickness to web 
t b 

spacing — , however, is varied from one wing to the other. 
b w 

The solutions are initially in agreement, but when the webs absorb 
a sufficient amount of heat to decrease appre siably the temperature 
gradient between the skin and webs, the two solutions diverge. The 
torsional stiffness, as obtained from the present analysis, reaches a 
minimum. On the other hand, the torsional stiffness reduction obtained 
from the solution permitting no heat transfer to the webs continues to 
decrease and approaches an asymptotic value as the skin temperature 
approaches the adiabatic wall temperature. 
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In evaluating the constant heat- transfer coefficient, the value of 
h at the midchord (d = c/2) was chosen as reference and evaluated by 
the use of equation (6) for turbulent flow. The following values cor- 
responding to an altitude of 50^000 feet were used in the calculations: 
The velocity of sound was taken to be 971 ft/sec; the kinematic viscosity 

was assumed to be 8.18 x 10-^ sq ft/sec; k a = J.l6 x 10"^ Btu/ft-sec-°F; 
l/d 

and Np r = 0-9* The value of c m p for stainless steel was taken as 
59.1 Btu/ft5-°F. 

Effect of Variations in Heat-Transfer Coefficient 

Effect of turbulent and laminar boundary- layer flow .- In the pre- 
vious section the heat-transfer coefficient was assumed constant along 
the chord. Heat-transfer coefficients under actual flight conditions, 
however, vary along the chord as can be seen from reference 2. In order 
to study the effects of a variation in heat-transfer coefficient on the 
reduction in torsional stiffness, h was assumed to vary continuously 
across the chord in the manner given by equations (7) and (9) for tur- 
bulent and laminar boundary- layer flow. In equations (7) and (9) h 
was expressed in terms of the heat- transfer coefficient at the midchord 
in order to obtain a reference for comparison purposes. Inasmuch as the 
heat-transfer coefficients are nonsymmetrical about the midchord, the 
temperature distributions are no longer symmetric and consequently both 
a translation and a rotation of the cross section is possible. Further- 
more, with the heat-transfer coefficient a function of chordwise posi- 
tion, an analytic solution for the torsional-stiffness parameter $ is 
rather difficult to obtain. For the present analysis a numerical inte- 
gration procedure described in appendix C was employed. This procedure 
was used to obtain the reduction in torsional stiffness for two specific 
wing configurations subjected to variable heat- transfer coefficients 
across the chord. The results are presented in figures 9 and 10. 

In figure 9 the reduction in torsional stiffness of a solid double- 
wedge wing in turbulent and laminar boundary- layer flow is shown. For 
the sake of comparison the curve of torsional- stiffness reduction for 
a constant value of h is also presented. In the calculations of these 
results 20 chordwise stations were used. This number of stations was 
necessary to give accuracy within 1 percent when the calculations were 
compared with the analytical solution presented for a constant value of 
h in reference 1 . 

As shown in figure 9, the variation of heat-transfer coefficient 
does not substantially influence the magnitude of the torsional- stiffness 
reduction. The percent difference between $ calculated for a variable 
and for a constant heat-transfer coefficient is 7 percent for laminar 
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flow and only 3 percent for turbulent flow. Thus, the stiffness reduc- 
tions due to varying the heat- transfer coefficient are but a small per- 
centage of the stiffness loss resulting from the mass distribution. 

In figure 10 results are presented for the reduction in torsional 
stiffness of a hollow double-wedge airfoil in turbulent and laminar 
boundary- layer flow. For the hollow double wedge of constant skin thick- 
ness the reduction in torsional stiffness is due solely to chordwise 
variations of h. In the calculations of these results 40 chordwise 
stations were used. This number of stations was necessary to give accu- 
racy within 1 percent when the calculations were compared with an exact 
solution derived in reference 1 for turbulent flow over a hollow double- 
wedge cross section. Note that the maximum value of $ calculated for 
the hollow double wedge with laminar boundary- layer flow is only 38 per- 
cent of that calculated for a solid double wedge. 

Effect of a step transition in heat-transfer coefficient .- In order 
to investigate the effect on stiffness reduction of a transition in heat- 
transfer coefficient along the chord a simplified approach to the problem 
has been made. The heat-transfer coefficient is assumed to vary as a 
step function along the chord. Because of the difficulty in obtaining 
an analytic solution, numerical integration procedures described in 
appendix C were employed. Calculated results are presented for step 
transitions at the midchord and quarterchord of a solid double-wedge 
section. From a comparison of the results obtained by the numerical 
procedure and those obtained by an analytic solution derived for a step 
transition at the midchord, it was found that 20 chordwise stations were 
sufficient to give agreement within 2 percent. 

For a step transition at the midchord, tl 
any point along the chord can be expressed by 

-AC 

0 = 1 - e 1+ ^ 2x / C ' 1 
-A 

e . 1 - e 1 -^" 

where 

C = = Constant 

*He 

? _ 2h te T 

Cfflpt mftY 


he average temperature at 


(-c/2 5 x < 0) 
(0 < x % c/2) 


( 22 ) 


(23) 


( 24 ) 
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For a step transition at the quarterchord, the average temperature at 
any point along the chord is written as 


-AC 


= 1 - e 1+(2x / c) (-c/2 £ x < -c/4) 


e = 1 - e 


0 = 1 - e 


-A 

l+(2x/c) 

-A 

l-(2x/c) 


(-c/4 < x < 0) 


(0 < x S c/2) 


? (25) 


The thermal stresses and torsional stiffness reductions are eval- 
uated from the general derivation given in appendix C. 


In figure 11 $ is plotted as a function of A for a step transi- 
tion at the midchord for several different ratios of heat-transfer coef- 
ficient. For all the cases presented, the heat-transfer coefficient in 
the rearward portion remains the same and A is evaluated using this 
heat-transfer coefficient. The maximum stiffness reduction occurs in 
the limiting case where the heat- transfer coefficient is of constant 
magnitude across the entire chord (C = l) . For progressive reductions 
in the leading-edge heat-transfer coefficient, the reduction in tor- 
sional stiffness decreases. 


The effect of varying the point of transition along the chord is 
shown in figure 12 for a ratio of h le /h te = l/5- The transition point 

is moved from the leading edge, to the quarterchord, to the midchord. 
For the three transition points considered, the smallest reduction in 
torsional stiffness occurs with the step transition at the quarterchord 
while the leading-edge position results in the greatest reduction. 


CONCLUDING REMARKS 


The distribution of the heat-absorbing mass of an airfoil cross 
section is the primary factor in determining the reduction in torsional 
stiffness induced by aerodynamic heating. 

The percent reduction in torsional stiffness in a hollow double- 
wedge section, which is due entirely to the chordwise variation in heat- 
transfer coefficient caused by laminar and turbulent flow, is much 
smaller in magnitude than that calculated for a solid double wedge . 
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Calculations based, on the solid double wedge indicate that the variation 
in heat-transfer coefficient due to turbulent and laminar boundary- layer 
flow had very little effect on the magnitude of the stiffness reduction. 

From the results presented for the modified solid double wedge with 
a constant heat-transfer coefficient, it can be concluded that the flat- 
tening of the cross section reduces the maximum stiffness losses suffered 
by the unmodified wedge. On the other hand, slight degrees of blunting 
increase the torsional— stiffness reduction. It has been shown that, for 
both the blunted and the blunted and flattened wedge, ratios of tip thick- 
ness to midchord thickness up to approximately 0.20 exhibit an increase 
in torsional- stiffness reduction over that for the unblunted wedge. The 
percent variation in torsional-stiffness loss due to blunting seems to be 
independent of the magnitude of flattening of the cross section. 

For the case of a constant value of heat-transfer coefficient the 
loss of torsional stiffness occurring in multiweb wings of light con- 
struction was substantially less than that of corresponding solid double- 
wedge wings . 

The effect of a chordwise step transition in heat- transfer coeffi- 
cient on the torsional- stiffness reduction of solid wings was found to 
depend on the magnitude of the resultant jump and upon the point of 
transition. In all the cases studied in this investigation, the influ- 
ence of the step transition was to decrease the torsional- stiffness 
losses over those obtained with a uniform value of heat— transfer 
coefficient . 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., October lb, 1S>58* 
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APPENDIX A 

REDUCTION IN TORSIONAL STIFFNESS OF A MODIFIED DOUBLE WEDGE 


The modified double-wedge section to be analyzed is shown in fig- 
ure 1. For the purpose of this analysis, the assumptions mentioned in 
the section entitled "Basic Theory" are used. In addition, the heat- 
transfer coefficient is assumed constant along the chord. From equa- 
tion (4) then, the temperature distribution which Is symmetrical about 
the midchord can be written as 


e = 1 - e 


-A 


(0 x ^ a) 


e = 1 - e 


-Ab 

a+b-x 


( aSxS f) 


> (Al) 


With symmetry of temperature distribution about the midchord, the span- 
wise stress can be expressed by equation (ll) where / 0x dA = 0 as 

J A 


0 

Ea AT aw 



0 


(A 2) 


Substituting equation (Al) into equation (A2) and integrating yields 


-Xb 

<J _ g a+b-x ^max 


Ea AT~ 


Atb 


-Xb 

-e'^b 2 - Xb 2 + 2b a) + e L (l 2 - Xbl.) + X 2 b 2 F 0 ^a $ x % | j 


> (AJ) 


o _ ^max 
Ea AT aw A^b 


-A 


-Ab 

(-Xb 2 + L 2 ) + e L (l 2 - XbL^ + X 2 b 2 F Q 


(0 < x S a) 
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where F 0 = Ei(-A) - Ei^j, L = a + b - |, and Ei(-A) is the expo- 
nential integral function defined by 


Ei(-A) = - f - — dp (A4) 

'-'A P 


If the expressions for stress in equation (A;) are substituted into 
equation (l) and the assumption is made that J Q = ^Ixx> reduction 

in torsional stiffness becomes 


QJeff 

GJ 0 
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Lt + ^ a "^max ) 


r 

Jo 




where 


L 


c / 2 0i (a + b - x) 
*** ^max T 


cbc 


= |^max 5b - Lt z e 5 + ^ at max 5 ^ 


(A5) 


(A6) 


Integrating and simplifying the above equation results In the following 
expression: 


QJeff 

GJ 0 



where 



12 



v. 
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APPENDIX B 


REDUCTION IN TORSIONAL STIFFNESS OF A MULTIWEB WEDGE 

The multiweb wing to be analyzed is shown in figure 6. The assump- 
tions previously mentioned in the section entitled "Basic Theory" will 
be employed with the heat-transfer coefficient considered constant in 
the chordwise direction. 

The temperature distribution of each skin of the multiweb is assumed 
to vary as that of a flat plate; hence. 


-A 


e s = 1 - 


^tg/tjnax 


(Bl) 


The temperature distribution of the web material is assumed to vary 
as if the webs completely filled the interior of the multiweb and formed 
a solid double wedge. Thus, 


®w ~ 1 



(B2) 


The average temperature in a typical bay of the multiweb can be written 
as 


2^w"^s®s * ttryfBv 

2b w t s + tt” 


(B.J) 


In order to evaluate the spanwise stress at any point along the 
chord, the total web material is considered to vary continuously along 
the chord in proportion to the wing thickness. The expression for 9 
given in equation (B3) is then substituted into equation (ll) for span- 
wise stress. Because of symmetry of temperature distribution about the 


midchord, 


0x dA 


1 A 


in equation (ll) is equal to zero. 


Since 
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dA = ct c 


1 + d.0^ , by letting ^ | the expression for span- 


wise stress becomes 


en + 2i ^jdi 


Ea AT 


aw 


'o 1 + 2l 5j dt 


- e 


(b4) 


where 


Thus, 


Aw tjoaxW 

As b\t s 


Ea AT av 


1 + 


*w 


fo ( 6s + 2|9v ^) d| ' 6 


(B5) 


Substitution of equation (B5) into equation (l) yields, 



The initial torsional stiffness of the multiweb GJ 0 is derived 
by assuming that the webs are rigid in shear. Thus, 


Jo 



2 

max 


t s 


(B7) 


By integrating and simplifying equation (B6) , the reduction in torsional 
stiffness can be expressed as 
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where 


_ ~1 J ge^s/^inax/ ^ _ 


A s + A w 


— e _A A 5 + JX 2 + 4A + (1 - A) 


Aft + At 


A 2 ^ fJa 2 + 8A + 10 


Aa + At. 


j ->co _p 

- — d.p = Ei(-A) 
A ^ 
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APPENDIX C 

EFFECT OF A VARIABLE HEAT-TRANSFER COEFFICIENT ON 
REDUCTION OF TORSIONAL STIFFNESS 


In appendixes A and B the heat-transfer coefficient was assumed to be 
constant and the expressions for spanvise stress and torsional-stiffness 
reduction could he easily expressed in a closed form. However, when t e 
heat-transfer coefficient is varied along the chord, a numerical integra- 
tion procedure proves to be more expedient. 

Consider any arbitrary symmetric cross sedtion with a variable 
heat-transfer coefficient across the chord. If the cross section is 
first divided into N number of equally spaced elements, the expression 
for spanwise stress in equation (ll) can be written as 





I x i 2 n 

i=l 


e i 


where 

0^ temperature of each element 

t± thickness of each element 

Xi distance of each element from midchord 

Substitution in equation (l) yields 


(Cl) 


OJeff 

GJ 0 



N 


^ <*i x i 2t i 


i=l 


(C2) 



2k 


where 8 is the width of an element. Let 


Qi = — = — 
Ea AT aw 

Then equation (C2) can be expressed as 


(C3) 


GJ. 


eff „ Ea 


GJ r 


1 AT 


awl 


L Wx / 


where 


N 


* ' r I xl2tiQ i 


1=1 
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Figure 7.- Effect of area distribution of multiweb on torsional - 

stiffness reduction. 



Figure 8.- Comparison of multiweb solutions. Thickness ratio, 3 percent 
Mo = 0.75 - M-f = 4.0; altitude, 50,000 feet; — = 0.35* 
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- Effect of a variable heat-transfer coefficient on torsional- 
stiffness reduction of a solid double-wedge wing. 



Figure 10.- Effect of a variable heat-transfer coefficient on torsional 
stiffness reduction of a hollow double-wedge wing. 




Figure 11.- Effect of transition at midchord on torsional-stiffness 
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Figure 12.- Effect of varying transition-point on torsional-stiffness 
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